PACIAE 2.0: An updated parton and hadron cascade model (program) for 
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We have updated the parton and hadron cascade model PACIAE for the relativistic nuclear 
collisions, from based on JETSET 6.4 and PYTHIA 5.7 to based on PYTHIA 6.4, and renamed as 
PACIAE 2.0. The main physics concerning the stages of the parton initiation, parton rescattering, 
hadronization, and hadron rescattering were discussed. The structures of the programs were briefly 
explained. In addition, some calculated examples were compared with the experimental data. It 
turns out that this model (program) works well. 
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Title of program: PACIAE version 2.0 
Catalogue number: ADBS 

Program obtained from: CPC Program Library, Queen's University, Belfast, N. Ireland; 
or from yanyl@ciae.ac.cn, zhoudm@phy.ccnu.edu.cn, sabh@ciae.ac.cn 

Computer for which the program is designed and others on which it has been tested: DELL Studio XPS and 
others with a FORTRAN 77 or GFORTRAN compiler 

Computer: DELL Studio XPS 

Programming languages: FORTRAN 77 

Memory required to execute with typical data: ~1G words 

No. of bits in a word: 64 

Peripherals used: terminal for input, terminal or printer for output 
No. of lines in distributed program, including test data, etc.: ~200000 
Distribution format: tar.gz 

Nature of physical problem: The Monte Carlo simulation of hadron transport (cascade) model is successful in 
studying the observables at final state in the relativistic nuclear collisions. However the high pt suppression, 
the jet quenching (energy loss), and the eccentricity scaling of vi etc., observed in high energy nuclear collisions, 
indicates the important effect of the initial partonic state on the final hadronic state. Therefore the better parton 
and hadron transport (cascade) models for the relativistic nuclear collisions are highly required. 

Method of solution: The parton and hadron cascade model PACIAE is originally based on the JETSET 7.4 
and PYTHIA 5.7. The PYTHIA model has been updated to PYTHIA 6.4 with the additions of new physics, 
the improvements in existing physics, and the embedding of the JETSET model etc.. Therefore we update the 
PACIAE model to the new version of PACIAE 2.0 based on the PYTHIA 6.4 in this paper. In addition, some 
improvements in physics have been introduced in this new version. 

Restrictions: Depend on the problem studied. 

Typical running time: Depend on the type of collision and energy. Examples running on the DELL Studio 
XPS are follows: 

• Running 1000 events for inelastic pp collisions at ^=200 GeV by program PACIAE 2.0a to reproduce 
PHOBOS data of rapidity density at mid-rapidity, dN c h / ^2/=2.251q'|q spends w 3 minutes. 

• Running 0-6% most central Au+Au collision at ,/i^"=200 GeV by program PACIAE 2.0b and PACIAE 
2.0c to reproduce PHOBOS data of charged multiplicity of 5060 0] spends «13 seconds/event and ~265 
seconds/event, respectively. 

PACS: 25.75.Dw, 24.10.Lx 

Keywords: relativistic nuclear collision; transport (cascade) model; hadron; parton; parton rescattering; 
hadronization; hadron rescattering. 
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LONG WRITE-UP 



I. INTRODUCTION 

The hadron transport (cascade) model is successful in describing the relativistic nuclear collisions. 
However, many new phenomena observed in the relativistic nucleus-nucleus (including proton-nucleus) 
collisions at RHIC energies, such as highp T suppression 0,0], jet quenching (energy loss) @, and elliptic 
flow eccentricity scaling jfij etc., strongly indicate the important effect of the partonic initial state on the 
hadronic final state. The better parton and hadron transport (cascade) model is urgently required. 

The first parton and hadron cascade model encouraged copious publications of the similar models, 
such as AMPT (string melting) [3, PACIAE H, BAMPS (parton cascade only) [13], PHSD [TO], and 
MARTINI [12;] etc.. Among them the PACIAE model, a parton and hadron cascade model for relativistic 
nuclear collision, is based on PYTHIA 5.7 and JETSET 7.4 [Hj]. We have employed this PACIAE model 
successfully studying [14j: 

• limiting fragmentation phenomenon in elementary and heavy-ion collisions, 

• fluctuations and correlations in particle production, 

• direct photon production in pp and Au+Au collisions at RHIC energies, 

• property of QCD matter at RHIC and LHC energies, 

• charged particle and strange particle productions in pp collisions at LHC energies, 

• charged particle elliptic flow in pp collisions at LHC energies. 

In this paper, the PACIAE model is updated to PACIAE 2.0 based on PYTHIA 6.4 with the additions 
of new physics and the improvements in existing physics. As the latest version PYTHIA 8.1 [l6[ "does 
not yet in every respect replaces the old code" and the physical differences between PYTHIA 8.1 and 
PYTHIA 6.4 are not important for the investigation of nucleus-nucleus collisions, so we left the connection 
with PYTHIA 8.1 in the future. 

PYTHIA 6.4 is a Monte Carlo event generator for relativistic hadron-hadron collisions in hadronic 
level. In this model a pp (hadron-hadron, hh) collision is decomposed into parton-parton collisions. A 
hard parton-parton collision is described by the lowest leading order perturbative QCD (LO-pQCD). 
The soft parton-parton collision, non-perturbative phenomenon, is considered empirically. Because the 
initial- and final-state QCD radiations and multiparton interactions are considered in the parton-parton 
scattering, the consequence of a hh collision is a parton multijet configuration composed of di-quarks 
(anti-diquarks), quarks (anti-quarks), and gluons, besides a few hadronic remnants. This parton multijet 
configuration is followed by the string construction and fragmentation (hadronization). Therefore one 
obtains a hadronic final state for a hh (pp) collision. If one switches off the string fragmentation and 
breaks up the di-quarks (anti-diquarks) to quarks (anti-quarks), one then obtains a state composed of 
the quarks, anti-quarks, and the gluons. This is the key point in the creation of the parton and hadron 
cascade model PACIAE for the relativistic nuclear collisions. 

The PACIAE model composes of four stages of the parton initiation, parton rescattering, hadronization, 
and hadron rescattering. PACIAE 2.0 has three versions: PACIAE 2.0a describing the relativistic pp 
collision (pp, or e + e~, denoted as elementary collision later) and PACIAE 2.0b as well as PACIAE 2.0c 
describing the relativistic nucleus-nucleus (A+B, including p+A) collisions. The hybrid system of units, 
based on the mass in GeV, length in fm, and time in fm with c=l and hc=1.2A GeV-fm, is used in 
PACIAE 2.0. 

The rest of this paper is organized as follows: Section II discussed the main physics in the PACIAE 
model except the ones in PYTHIA: 

• nucleon initiation in the position and momentum phase spaces, 

• particle propagation (cascade) and the collision criteria, 

• cross sections, 

• determination of the scattered particles, 

• diquark break-up, 
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• reduction of the strange (heavy) quark suppression, 

• deexcitation of the energetic quark (antiquark) in the coalescence model. 

We explain the structures of programs in section III. The examples are calculated and compared with 
the experimental data in section IV. We give the users' guide in the appendix. 

II. PHYSICS CONCERNED 

The PYTHIA 6.4 model involves abundant particle and nuclear physics contents which have been 
discussed in fl5l j in the detail. Here we just introduce the extra physics concerned with the nuclear 
initiation, the parton rescattering, and the hadron rescattering. 

A. Nucleon initiation in the position and momentum phase spaces 

1. Impact parameter and its sampling 

In the experiments the centrality bins are defined by the cuts in the particle multiplicity distribution 
and indicated by the percentages, g, in the geometric (total) cross section. However, the centrality is 
conveniently defined theoretically by the impact parameter b. A mapping relation between g and b 

b = ^b max , b max = Ra + Rb (1) 

was introduced according to the definition of geometric cross section 

As the experiments of nucleus-nucleus collisions are always performed in a given centrality bin the 
relevant theoretical calculations should span the corresponding impact parameter interval from bd to b u . 
We sample impact parameter bi randomly by 



h = ^x{bl-bl) + b% (2) 

and generate a nucleus-nucleus collision event. In the above equation £ refers to a random number. The 
observable (operator) averaged over the events generated by different bi (i=l,2,...) can compare with the 
corresponding experimental data. 

One may also use the systematic sampling method: 

fb 2 rb n =b u ,2 _ u2 

bdb= bdb=...= bdb= " d , (3) 

b =b d Jbl Jb n -l 2rl 



b j = \li b -^-^+ b l j = Q,l,2,...,n(neven) (4) 



n 



The set of bj(j = 1, 3, 5, n— 1) is then the sampled spectral values. We generate events by these spectral 
values repeatedly and compare the event averaged observable with the corresponding experimental data. 

2. Geometric overlap zone and number of participant nucleons 

For a nucleus- nucleus collision A+B at a given impact parameter b, a geometric overlap zone is formed 
between colliding nuclei in the position phase space. The nucleons located inside this zone are the 
participant nucleons, otherwise the spectator nucleons. The geometric number of participant nucleons is 
calculated [ijj by 

N part (b) =N p A art (b) + N* rt (b) , (5) 
N* rt (b) = PA J dV6(R A - (x 2 + (b- yf + z 2 )^ 2 )9(R B - (x 2 + y 2 ) 1 ' 2 ), (6) 

N p B art (b) =p B f dV9(R B (x 2 +y 2 + z 2 ) l ' 2 )6{R A - (x 2 + (b - y) 2 ) 1 / 2 ), (7) 
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where 

®^ | 1 otherwise ^ 

Ra and pa (Rb and ps) are the radius and nuclear density of the nucleus A (B), respectively. The nuclear 
density of nucleus A is normalized to the atomic number A (we denote the nucleus and its atomic number 
by A uniquely). The STAR and PHOBOS collaborations have reported that their analysed Glauber model 
number of participant nucleons is nearly reaction energy independent for a given centrality A+B collisions 
within error bars [l], [TH- Therefore we assume normal nuclear density of pa = Pb = Po— 0.16 fm~ 3 here. 
The geometric number of spectator nucleons is then calculated by 

spec parti spec part' \ ) 



3. Nucleon initiation in the position and momentum phase spaces 



The geometric participant nucleons in each of the colliding nuclei are distributed randomly inside the 
geometric overlap zone. And the geometric spectator nucleons are distributed according to the Woods- 
Saxon distribution and A-k uniform distribution and requiring to be outside the geometric overlap zone. 
These distrbutions are 

p{r) = pa {\ + C ^ r —J^y\ (10) 

/('.*) = («) 

where R = ro^4 1 ^ 3 +0.54 fm is the radius of nucleus A, the stander root-mean-square-radius parameter 
ro ~1.15 fm, and d ~ 0.54 fm is the length of diffusion tail. 

For the case of collider the x and y components of momentum of each nucleons in colliding nucleus 
are assumed to be zero, provided the beam is orientated towards the z direction. The beam momentum 
is given to p z for each nucleons properly. This means the effect of Fermi motion is neglected in the 
relativistic nuclear collisions. 



B. Particle propagation (cascade) and the collision criteria 

The particles (partons in the parton rescattering stage or hadrons in the hadron rescattering stage) 
are assumed to travel in the straight trajectories 

n = R l + p l {t~T l ) (12) 

along the momentum direction during two consecutive collisions. In the above equation (i?^, Ti) and Pi 
stand for the four-position and velocity of the i-th particle at the moment of last collision, respectively. 

The particles interact with each other by means of total cross section. In order to perform the collisions 
in a consistent way one needs to introduce a time ordering procedure. We calculate the time between 
the last and the next possible collisions for all particle pairs and then choose the minimum collision time 
pair to perform the particle-particle collision (ljjj . A particle-particle collision (particle i bombards with 
j, for instance) is defined to occur provided the minimum approaching distance (D) between particles i 
and j, in their cms system, satisfies [l9j 



D < y/<T tot /n, (13) 

where a tot is the total cross section in fm 2 . We calculate D starting from the trajectories of particle i 
and j in their cms system (variables in the cms system indicated with superscript *) 

r* =R* + 0*(t*-T*) , r*=R*+0*(t*-T*), (14) 
The squared relative distance between i and j reads 

d 2 (t*) = (AR*) 2 + 2AR*Aj3*[t* ^(/3*T* - 0TJT)] + (A/3*) 2 [i* —0ST7 ~ P*T*)} 2 , (15) 

A/3* J J ' 
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*« = /A ;.„;? 3 1 1 - (is) 



where 

AR* = R*-R*, A/3* =p*-0*. (16) 

From the requirement of 

ffi = (17) 

at* v ; 

one obtains the collision time at the moment of i bombarding with j 

-Ai?*A/3* + A/3*(/3*T* - f3*T x 
(A(3*) 2 

Inserting Eq. (|18[) into Eq. (|15p . the corresponding minimum approaching distance is obtained 

D = [(AR*) 2 + AR*Af3*(t*A -(/3*T* - PIT*))] 1 ' 2 . (19) 

As the particle collision must be causal, the above calculated collision time should satisfy 

t*. >max(:77,T*). (20) 

Eqs. (fT3")) and (|20p are the criteria for a particle-particle collision to occur. 

This collision time has to boost back from the cms of colliding pair to the cms of nucleus-nucleus 
collision according to the inverse Lorentz transformation [20L l2lj 

P 2 

t =j(t* + f*-p), 

3 = Pj±Pj ( 21 ) 
' E t + Ej ' 
1 

where /3 refers to the velocity of cms of colliding pair in the frame of the cms of nucleus-nucleus collision. 
Ordering the collision times calculated for all collision pairs, the collision (time) list is obtained. 
The corresponding Lorentz transformation reads [2(| Hl[ 

r* = r + p[ — 2 jt], 

t*=~f(t-r-0). (22) 

Replacing four position (r, i) with four momentum (p, E) in Eq. (|22p and (1211) one obtains, respectively, 
the Lorentz transformation and inverse Lorentz transformation for the four momentum. 



C. Cross sections 

1. Hadron-hadron cross sections 

The isospin averaged parametrization formulas are employed for the hh cross sections 0, [23[ ■ Taking 
the 7r + p — > k + Y (Y refers to A and E) reactions as example, the cross sections are 
n+p -> E+fc+: 

' Ai(Vs- 1.683), 1.683 GeV < Vs < 1.934 GeV, 
A x = 0.7 mb/0.218 GeV, 

A 2 exp(-A 3y /s), 1.934 GeV < y/s < 3.0 GeV, 
a(y/s) = { , (23) 

* A 2 = 60.26 mb, A 3 = 2.31 GeV , V ' 

A A exp(-A 5 y/s), 3.0 GeV < y/s, 

A 4 = 0.36mb, A 5 = 0.605 GeV' 1 , 



G 



7T p 



Ak°: 



7T p 



' Ai(y/s- 1.613), 1.613 GeV < Vs < 1-684 GeV, 
A x = 0.9 mb/0.091 GeV, 

A 2 exp(-A 3 Vs), 1-684 GeV < Vs < 2.1 GeV, 

A 2 = 436.3 mb, A 3 = 4.154 GeV _1 , 

A 4 exp(-A 5V / s), 2.1 GeV < Vs, 

A A = 0.314 mb, A 5 = 0.301 GeV" 1 , 



' Ax{yfs - 1.689), 1.689 GeV < Vs < 1.722 GeV, 

Ax = 10.6 mb/ GeV, 
A 2 exp{-A 3 y/H), 1.722 GeV < Vs < 3.0 GeV, 

A 2 = 13.7mb, A 3 = 1.92 GeV" 1 , 
Ax exp(-A 5 Vs), 3.0 GeV < Vs, 

A 4 = 0.188 mb, A 5 = 0.611 GeV -1 , 



A 6 (constant), 1.691 GeV < Vs < 1.9 GeV, 
A 4 exp(-A 5 Vs), 1.9 GeV < Vs, 

A 4 = 309.06 mb, A 5 = 3.77 GeV" 1 . 



(24) 



(25) 



(26) 



An assumed constant total cross sections of = 40 mb (RHIC energy region) or 76 mb (LHC energy 
region) [2l|, — 25 mb, — 21 mb, and er^JJ = 10 mb as well as the inelastic to total cross section 
ratio of 0.85 are provided as another option. We assume 



pp nn pn _ An AA 

'tot — "tot — "tot ~ "tot — "tot ■ 



(27) 



The above hh cross sections are used in the hadron rescattering stage and in the parton initiation stage 
for the nucleus-nucleus collisions as well. 
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2. Parton-parton cross sections 



We consider only the 2— > 2 parton-parton (re)scattering at the moment. The nine parton-parton 
differential cross sections calculated by means of LO-pQCD 0, HH are: 



— (aft -> cd:s,i) = — ^ M(a6 -> cd) ' 
tit s z 



\M{q l q l -> qiqi)\' 



\M{qiqj -> 

4 s 2 + u 2 s 2 + £ 2 , 
1 9 ( ~ + u 2 ' 

\M{qiqj -> q t q 3 )\ 2 - ■ 



4s 2 +w 2 _ 8 s 2 



9 t2 



27fi£ ~ 9P 



4 s 2 + u 2 



\M(qiqi -> q iqi )Y = -( 



4 / s 2 + u 2 i 2 + u z 



9 V f2 

ITTF, M2 32 u 2 +£ 



\M{gg -> 9i $)| 2 
|M((z l5 -^<z. i3 )| 2 



4 g + M 2 
9 s 2 ^ 

27 §t ^ 9 t 2 
.u 2 + i 2 ^ 32 s 
27 ^| 3 s 2 " ~27i 
1 u 2 + i 2 3 u 2 + t 2 ^ Is 
~ 6 7i£ 8 J 2 ~ _ 3f 
J' 

^ I 2 



9 us 



i 1 



ut us si 



9 s 2 



\M(gg -> gg)\ 2 = -(3 - ^ - it - tj) « -57 



< 2 u 2 



2# 



(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
(34) 
(35) 
(36) 
(37) 



Here the parton is represented as a classical point-like particle and specified by its internal degree of 
freedoms: flavor, rest (current) mass, momentum, and position. The spin and color degrees of freedoms 
are not explicitly included [7J. Therefore the amplitudes \M\ 2 are calculated by averaging over spin and 
color degrees of freedoms. The parton energy is determined by E 2 = p 2 + m 2 + M 2 , where M refers 
to virtual mass (M 2 =0 for on mass shell parton, M 2 <0 for space-like parton, and M 2 >0 for time-like 
parton) 0]. a s refers to the effective strong coupling constant, s, i, and u stand for the Mandelstam 
variables. In each of the above equations the second form is obtained by keeping only the leading divergent 
terms [26j with the assumption of the parton current mass is negligible relative to the cms energy. This 
assumption leads to 



8 + 1 + u = Q. 



(38) 



The subprocesses of Eqs. (|32|) . (|34|) . and (f35)) are inelastic (re)scatterings, otherwise elastic. Subprocess 
of Eq. (f3"2"|) is always negligible. 

If the thermal direct photon is interested following two subprocesses have to be included: 



16 a 



\M(qg->qj)\ 2 

where a=l/137 stands for the fine-structure constant, e q refers to the charge of quark q. 



1 a 
3 a 



. .. 9 , i s . la., 

o— 4 - + t ~-o— e h 



(39) 
(40) 



D. Determination of the scattered particles 



The particle-particle (re)scattering is always dealt in the cms of colliding pair. The superscript * is 
neglected in the rest of this section, except mentioned specially. 
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1. Hadron-hadron collisions 



For an inelastic hadron (re) scattering 1 + 2 — > 3 + 4 (mi ^ 7713, mi 7^ m.4), the four momentum of the 
scattered particle reads 



_. _> _ r - (m 3 + m 4 ) 2 ][s - (m 3 - m A y 



\PZ\ Z =\P4\ 



4.s 



£3 = i± (mf-m|) ! 
s — (m 3 — mf) 

£ 4 = ^— F • 



They are deduced from four momentum conservation |20t |21| . In the above equation s refers to the cms 
energy squared. 

The momentum orientation (polar angle and azimuthal angle <fi) of scattered particle is decided by 
the assumption that the momentum transfer squared, t, is distributed as [19|, |27[ 

da „ , , , 

— - exp Bt, (42) 

= 3.65(vs — mi — 7712) ~ 3.65\/s, (43) 



o 



a 6 



* = _A«4 (44) 

A = min(10.3, (1.12 < p T >) -2 ), (45) 

where < px > stands for the mean transverse momentum assumed to be 0.5 GeV/c because it is observed 
in the experiment that the value of < pt > may be not so sensitive to the reaction energy and centrality 
[28l ]. In the above equations the second form was obtained provided the particle mass is negligible for 
high cms energy. On the other hand, t is related to 9 S spanned between p\ and JJ3 [2l| by 

t = m\ - 2eie 3 + 2pi • p s + m 2 3 , (46) 
Unax = m\ - 2ex€ 3 + 2|^i||p 3 | + m\, (47) 
tmin = mf - 2eie 3 - 2|pi||p 3 | + m\, (48) 

where t max and t m i n are corresponding to cos0 s =l and -1, respectively. Therefore we first sample a t' 
value within [i TO , n , t max ] according to exponential distribution Eq. (|42|) 

t' = - ln[eexp(Br; max ) + (1 - exp(Bt min )]. (49) 
t> 

The corresponding cos 6 S is then calculated by Eq. (|4"6")l 

co S g s =°- 5(t/ - m j- wi)+ei£3 . (50) 

|Pl||P3| 

The azimuthal angle <p s is sampled isotropically in 2n. 

It should be mentioned that the orientation of scattered particle 3, for instance, is relative to scattering 
particle 1, so it should rotate back to the particle 1 and 2 cms system, where particle 1 is described by 
\pi\, 0i, and ef>i. This rotation reads 



'AS J: 



p 3x =P3 X [cos <pi (cos 0i sin 9 S cos <p s + sin 0i cos S ) — sin cj>i sin S sin c 
p 3y =P3 y [sin 4>i (cos 0i sin0 s cost's + sin0i cos0 s ) + cos^i sin0 s sin0 s ], (51) 
p 3z ~P3z [cos 0i cos S — sin 0i sin S cos <j> 8 \ . 

At last the scattered particles are boosted back to the moving frame of nucleus-nucleus cms system. 
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It is impossible to take all inelastic channels into account, only following inelastic reactions 



TTN r 

NN 
nN ; 



7rA 

* iVA 
= kY 



irY r± kZ 
kN # ttY 
kY ^ ttS 
kN ^ fcS 
7rS ^ kVL~ 



TTiV 
7T7T 

nN 
ttY 
kN 
kY ^ nZ 
kN # fcS 



pN 
± kk 
kY 

7T? 



N N annihilation, NY 'annihilation , 

are considered. The rest is treated as elastic reactions in a sense. Even so, there are already ~600 inelastic 
channels involved. 

Taking nN (re)scattering as an example, there are channels of 



ttN -t ttA, nN pN, ttN -> kY, 



(52) 



therefore the relative probabilities of these channels are invoked to select one among them. The cross 
section of the reverse reactions are calculated by means of detailed balance. 

The final state of TV TV and NY annihilations are simply treated as five particle state through NN — > 
put — > 5n and NY — > k*u — > k + An, respectively. The cross section of NY annihilation is assumed to be 
1/5 of the NN annihilation. 

For elastic hadron (re)scattering because 



mi = 7713, % = 7714, 

Eqs. (|4*6")) -(|5T)1) reduce, respectively, to 



\Pl\ = \P2\ = \P3\ = \Pl\, Cl = £3, £2 = £4, 



t = 2\p 1 \ 2 (cose s -l), 



0. 



-4|pi| 2 , 



(53) 



(54) 
(55) 
(56) 



t' = - ln[£ + (1 - exp(Bt, 
is 



COS0* = 1 



2|PiP 



(57) 
(58) 



The rest is the same as that of inelastic (re)scattering. 



2. Parton-parton collisions 

If several partonic final states can be reached by a single partonic initial state the following relative 
probability is invoked to select one among them 

a(ab — > cd; s)/a a b(s) (59) 

where 

f° da 

a(ab — > cd; s) — / — -(ab — > cd; s,t)di (60) 
J-g dt 

gqfe(s) = p{ab -» cd; s). (61) 

c.d 



10 



In the inelastic parton-parton (re)scattering the flavor of scattered quark (antiquark) is decided ac- 
cording to 

u : d : s : c • • • = j u : j d : -y s : j c ■ ■ ■ (62) 

where 7„, 74, j a , j c , ■ ■ ■ are input parameters and their default values are 1:1:0.3:10"" ••• . The c quark 
is not included presently. 

The momentum transfer squared t is sampled according to 

1 r l do 



o(ab -> ca; s) di' 



Then the (re)scattering angle is calculated by [20J, |21| 

2si 2t i , . 

cos0 s = l + — - — -—— ^ sa i + _ = l + ^_^, (64) 

[s - (m a + m b )\[s - (m a - mb)\ s 2\Pa\ 

where the second form is obtained provided the parton mass is negligible relative to the cms energy. The 
azimuthal angle is sampled randomly in 2tt. 



E. Diquark break-up 



In order to obtain the initial partonic state for nuclear collision system one has to switch-off string 
fragmentation in PYTHIA, to break-up the strings, and to break-up the diquarks in the strings. The 
diquark break-up is performed in its rest frame according to two-body decay kinematics [13, Hlj . The 
energy and momentum modulus of broken quarks are calculated by Eq. (1411) with and 7714 identified 
as the mass of broken quarks and s the squared four momentum of diquark. The orientation of one of 
the broken quarks is sampled isotropically in An and the other is just orientated oppositely 

In the above, the three momentum of broken quark is calculated relative to three momentum of diquark, 
so it must be rotated back to the frame where diquark is described according to Eq. (|51[) . It is also need 
to boost back to the moving frame of diquark according to inverse Lorentz transformation Eq. (|2"Tj) . 

One of the broken quarks is assumed to have the diquark four position. The other is assumed to be 
created at the same time but is arranged randomly around the first one within 0.5 fm in each of the three 
position coordinates. 

The gluons are also needed to be broken-up in the hadronization by the Monte Carlo coalescence 
model. Because of zero mass the gluon is broken according to three momentum conservation. The energy 
discrepancy between the breaking gluon and the broken quarks is shared among quarks and antiquarks. 
The flavor of broken quark is decided by Eq. (IB"2l . 



F. Reduction of the strange (heavy) quark suppression 



In the LUND string fragmentation regime the qq pair with quark mass m and transverse momentum 
Pt may be created quantum mechanically at one point and then tunnel out to the classically allowed 
region [15]. This tunnelling probability is given by 

, 7rm 2 , . -kv% , , , 

exp( exp— ^ (65 

where the string tension k is assumed to be ~1 GeV/fm csO.2 GeV 2 [15]. This probability implies a 
suppression of strange (heavy) quark production ofit:cf:s:cwl:l: 0.3 : 10~ n . The charm and 
heavier quarks are not expected to be produced in the soft string fragmentation process if there is no 
charm and heavier quarks in the string originally. They are expected to be produced only in the hard 
process or as a part of the initial- and final-state QCD radiations. The higher px quark pair is expected 
to be created provided the string tension is large. However, in the PYTHIA 6.4 model there are model 
parameters of 

• parj(l) is the suppression of diquark- antidiquark pair production compared with quark-antiquark 
production, 
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• parj(2) is the suppression of s quark pair production compared with u or d pair production, 

• parj(3) is the extra suppression of strange diquark production compared with the normal suppression 
of strange quark, 

• parj(21) corresponds to the width a in the Gaussian p x and p y transverse momentum distributions 
for primary hadrons. 

They are able to be tuned to reduce the strange quark suppression and to change the width of its pr 
distribution. 

We introduced a mechanism of the increase of effective string tension and hence the reduction of strange 
quark suppression in [29| . In that paper we assumed that the effective string tension increases with the 
increasing of the number and hardening of gluons in the string by 

K e " = Ko(l-0-°, (66) 



£= "° '* ■ (67) 

In^ + E^M^) 

In above equations, kq is the string tension of the pure qq string assumed to be ~ 1 GeV/fm. The 
gluons in the multigluon string are ordered from 2 to n-1 because 1 and n refer to quark and antiquark 
at two ends of string, kxj is the transverse momentum of gluon j with > sq and hrmax is the largest 
transverse momentum among the gluons. The parameters a=3.5 GeV and so=0.8 GeV were determined 
by fitting the hh collision data. It should be mentioned that Eq. (|67p represents the deviation scale of 
the multigluon string from the pure string. 

If A denotes parj(2) (parj(l), parj(3)) then by Eq. (|65f the following relations between two strings 
with, respectively, effective string tension of n\ and ar e obtained 

\2 = \W T , (68) 



oa = oiFf 77 ) 1/a , (69) 

where subscripts identify the strings. 

Above mechanism is involved in the PACIAE 2.0 model. Therefore one can first tune the parameters 
of parj(l), (2), (3), and (21) to fit the strangeness production data in a given nuclear collision at a given 
energy. Then the resulted parj(l), (2), (3), (21) and effective string tension can be used to predict the 
strangeness production in the same reaction system at different energy even in the different reaction 
systems. 



G. Deexcitation of the energetic quark (antiquark) in coalescence model 

If the energy of a quark (anti-quark) is large enough (e.g. >2 GeV), one should introduce deexcitation 
mechanism for this quark (antiquark) in the coalescence model. This is similar to the iterative approach 
of quark-antiquark pair generation in the LUND string fragmentation regime [13l Il5l |16J. Taking an 
energetic initial quark go as an example, a new q\q\ pair may be created from vacuum. This quark- 
antiquark pair brings a fraction z of energy (momentum) of go and leaves remnant of go behind. If the 
energy of g remnant is still large enough, a new g 2 g"2 pair creates and above processes repeats until the 
energy is not enough to generate a new quark-antiquark pair. 

In the case of positive longitudinal momentum of mother quark p az in order to obtain four momentum 
of the first daughter pair gigi we introduce the forward light cone variable 

W a = E +p 0z . (70) 

The LUND string fragmentation function 



f(z) (x z- l {l - z) a exp(-/3m^/z) (71) 
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with default values of a = 0.3 and (3 = 0.58 or the Field-Feynman parametrization 

f(z) oc 1 - a + 3a(l - z) 2 (72) 

with default value of a = 0.77 is used to sample randomly a fraction variable z\ for the first daughter 
pair qiq±. Its forward light cone variable is then 



Wt=ziW Q . (73) 



From Eqs. ([70} and ([75]). we obtain 



= 2 W ~ Wf } (75) 

Similarly, for the negative longitudinal momentum of mother quark, one has to introduce the backward 
light cone variable 



w = E - poz (76) 



and has the solution of 



E > = 1^ + ir } ' (77) 

- 4< ff ' - ^ »■ < 78 > 

for the first daughter pair qiqi- 

The pix and pi y of qiqi pair are sampled according to the two dimensional Gaussian distribution (30j 

exp[-p|/a]dp|, < p% < p Tma x, (79) 

with default values of a=0.5 (GeV/c) -2 , PTmax—^ GeV/c. Then we obtain 

Pis =pircos(0), pij, = pirsin(0), (80) 

with sampled isotropically in 27r. 

The gigi pair is then broken up in its rest frame according to two-body kinematics mentioned above. 
The flavor of quark (antiquark) is sampled according to Eq. (|52^1 . It is also needed first to rotate back 
to the frame where q\q~\ pair is described according to Eq. (|51[) and then to boost back to the moving 
frame of q\q\ pair according to inverse Lorentz transformation Eq. (f2lj) . 



III. PROGRAM (MODEL) STRUCTURE 



PACIAE 2.0 has three versions: PACIAE 2.0a (its program packets are compressed to 20a. tar. gz) 
describing the relativistic elementary collisions (pp, pp, or e + e~) and PACIAE 2.0b (20b. tar. gz) as well 
as PACIAE 2.0c (20c. tar. gz) describing the relativistic nucleus-nucleus (A+B, including p+A) collisions. 



A. PACIAE 2.0a structure 



The PACIAE 2.0a model is for elementary collisions, it differs from PYTHIA 6.4 in the addition of 
the parton rescattering before hadronization and the hadron rescattering after hadronization. In order 
to create the initial partonic state we switch off the string fragmentation temporarily in the PACIAE 
2.0 model. Then we obtain a partonic initial state (parton list) after the strings are broken up and 
the diquarks (anti-diquarks) are split up randomly. This partonic matter proceeds parton rescattering. 
After parton rescattering the partonic matter is hadronized by the string fragmentation or Monte Carlo 
coalescence model. The hadronic rescattering is followed. Therefore the PACIAE 2.0a, as well as PA- 
CIAE 2.0b and PACIAE 2.0c, consists of the parton initiation, parton evolution (cascade, rescattering), 
hadronization, and hadron evolution (cascade, rescattering) four stages. 

PACIAE 2.0a program consists of paciae_20a.f, parcas_20a.f, sfm_20a.f, coales_20a.f, hadcas_20a.f, and 
p_20a.f: 
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1. paciae_20a.f plays both functions of an example of user program and parton initiation. The later 
is performed by PYTHIA with the string fragmentation switched-off and the diquarks broken-up 
randomly. We obtain a configuration of quarks, anti-quarks, gluons, and a few hadronic remnants. 
This is a partonic initial state (parton list) for an elementary collision. 

2. parcas_20a.f performs the parton rescattering based on the parton list above. The key ingredients 
of this cascade process were: 

• The partons are assumed traveling on the straight trajectories along their momentum direction. 

• The partons interact with each other by the 2 — > 2 LO pQCD cross section otot [HI HH| or 
its variety of keeping only leading divergent terms. With the collision times of all possible 
interaction pairs the collision time list is composed. 

• A parton-parton collision pair with least collision time is selected from the collision time list 
and executed according to the LO pQCD differential and total cross sections. 

• Update the parton list by removing the scattering partons and adding the scattered (generated) 
partons. 

• Update the collision time list by removing the collision pairs containing any one of the scatter- 
ing partons and adding the new collision pairs constructed by one of the scattered (generated) 
partons with another one in parton list. 

• A new parton-parton collision pair with least collision time is selected from the collision time 
list. 

• Repeat last four items until the collision time list is empty (partonic freeze-out). 

3. sfm_20a.f executes the hadronization of partonic matter according to the LUND string fragmentation 
after the parton rescattering and string reconstruction by calling p20a.f. 

4. coales_20a.f performs the hadronization of partonic matter according to the Monte Carlo coalescence 
model after the parton rescattering. The key ingredients of this coalescence model are: 

• The energetic partons are first deexcited. 

• The gluons are forcibly split into qq pair randomly. 

• A hadron table, composed of mesons and baryons, inputted to the program. The pseudoscalar 
and vector mesons made of u, d, s, and c quarks, the B + , B°, B*°, as well as T are considered 
as mesons. The baryons include the SU(4) multiplets of baryons made of u, d, s, and c quarks 
(except those with double c quarks) and A°. 

• Two partons coalesce a meson and three partons a baryon (antibaryon) in hadron table ac- 
cording to the quark flavor structure of the coalesced hadron and the flavors, positions, as well 
as momenta of the coalescing partons. 

• If the coalescing partons can form either a pseudoscalar meson or a vector meson (e. g. ud 
can be a ir + or a p + ) the principle of less discrepancy between the invariant mass of coalescing 
partons and the mass of coalesced hadron is invoked to select one from them. The same 
principle is used in the baryon production (such as p and A + are composed of uud). 

• The momentum conservation is required. 

• The phase space constraint 

16tt 2 a o a , h 3 

— g— Ar Ap = — (81) 

is introduced as an option. In the above equation h 3 /d is the volume occupied by a single 
hadron in the phase space, c?=4 refers to the spin and parity degeneracies of hadron, Ar and 
Ap stand for the position and momentum distances between coalescing partons, respectively. 

One selects sfm_20a.f or coales_20a.f to hadronize the partonic matter by switch parameter adjl(12). 

5. hadcas_20a.f executes hadron rescattering after hadronization. This hadron cascade proceeds in the 
same way as the one in the parton cascade mentioned above. But one should note: 

• The hadrons of n, k,p, n, p(uj), A, A, S, S, 0, J/^l and their antiparticles are considered here 
instead of parton there. 
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FIG. 1: Structure of PACIAE 2.0b (left) and PACIAE 2.0c (right) for the dynamical simulation of a nucleus- 
nucleus (A+B) collision. 



• Instead of LO pQCD parton-parton cross section is the hh total cross section. 

6. p20a.f is different from PYTHIA 6.4 in the addition of mechanism for the reduction of strange 
quark suppression. 

7. usux.dat is an example for input file. Each program packet (paciae _20a.f, parcas_20a.f, ...) is 
flexibly to be modified and/or replaced. A run can stop at the end of any packet for different 
purposes. For instance, one just sets the switch parameter adjl(40)=l if one takes interest in the 
partonic output before parton rescattering. 



TABLE I: Charged particle pseudorapidity densities at mid-rapidity (\r)\ < 1) for the INEL pp collisions having at 
least one charged particle in the same region (INEL> 0^|<i) at y/s=0.9, 2.36, and 7 TeV measured by ALICE [32l |. 
The first data error is statistical and the second systematic. The corresponding results calculated by PYTHIA 
6.4 and PACIAE 2.0a were given as well. 
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FIG. 2: (Color online) (a) charged particle pseudorapidity distribution in INEL pp collisions at -^5=900 GeV 
measured by ALICE [33|] and compared with the PACIAE 2.0a and PYTHIA 6.4 results, (b) transverse momentum 
distribution where data taken from [341. 



B. PACIAE 2.0b structure 



PACIAE 2.0b program consists of paciae_20b.f, parini_20b.f, parcas_20b.f, sfm_20b.f, colaes_20b.f, had- 
cas_20b.f, and p20b.f: 

1. paciae_20b.f is a user program. 

2. parini_20b.f administrates the generation of an event for an A+B collision: 

• Initiation in the position and momentum spaces for an A+B collision: 

— One colliding nucleus is set on the origin in the position space and the other on the position 
apart from the origin by an impact parameter b in the x axis. Both nucleus centers are 
assumed having y = z = 0. The colliding nucleus is assumed geometrically as a sphere 
with radius Raib) f m - 

— The geometric number of participant nucleons is first calculated. Then the participant 
nucleons are arranged randomly in the overlap region of the colliding nuclei. The rest 
nucleons in the colliding nucleus are randomly arranged in the corresponding sphere ac- 
cording to the Woods-Saxon distribution (for radius r) and 47r isotropic distribution (for 
orientation) and are required to be outside the overlap region. 

— The beam momentum is given to the nucleons in the colliding nucleus properly. 

— An initial nucleon list for an A+B collision is then obtained. 



TABLE II: Total charged multiplicity in 0-6% most central Au+Au collisions at ^/snn~0.2 GeV @| and the 
charged particle pseudorapidity density at mid-rapidity (|?7| <0.5) in 0-5% most central Pb+Pb collisions at 
y / sjviv=2.76 TeV [33l ] as well as the corresponding results calculated by PACIAE 2.0b, 2.0c, and 2.0c_c (the same 
as PACIAE 2.0c but hadronized by the Monte Carlo coalescence instead of the string fragmentation). 



Reaction 


VsiViv TeV 


Exp. 


PACIAE 2.0b PACIAE 2.0c PACIAE 2.0c_c 


Au+Au 
Pb+Pb 


0.2 
2.76 


5060+250 1 " 
1601+60* 


4940 4961 4746 
1554 1542 1540 



taken from [2J. * taken from [351 . 



• The nucleons in one colliding nucleus interact with nucleons in the other one when they travel- 
ing along their straight trajectories according to the NN total cross section. The corresponding 
collision time is calculated for all possible collision pairs and then the initial NN (hh) collision 
time list is composed. 

• A NN collision pair with least collision time is selected from the NN (hh) collision time list. 
This NN (hh) collision, if its cms energy is large enough [y/s >4.5 GeV, for instance), executes 
by PYTHIA with the string fragmentation switched-off and diquarks broken-up. A parton 
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FIG. 3: (Color online) (a) charged particle pseudorapidity distribution in 0-6% most central Au+Au collisions 
at y/s^=200 GeV measured by PHOBOS @] and compared with the PACIAE 2.0b results, (b) transverse 
momentum distribution where data taken from [3f|). 



configuration (parton list) for this NN (hh) collision is then obtained. If the ener gy is not 
large enough the NN (hh) collision is dealt with the usual two-body elastic collision |19| . 

• A parton cascade, the same as the one mentioned in last subsection, proceeds for current NN 
(hh) collision by parcas_20b.f. 

• sfm_20b.f or coales_20b.f executes the hadronization of partonic matter in current NN (hh) 
collision by the LUND string fragmentation regime (calling p20b.f) or Monte Carlo coalescence 
model according to the switch parameter adjl(12). sfm_20b.f and coales_20b.f are some what 
different from sfm_20a.f and coales_20a.f in PACIAE 2.0a, respectively. p20b.f is different from 
PYTHIA 6.4 in the addition of the mechanism for the reduction of strange quark suppression 
and of the treatment for the charge conservation in fragmentation of a string. 

• hadcas_20b.f performs hadron rescattering after hadronization of current NN (hh) collision. It 
is similar to the hadron rescattering in PACIAE 2.0a. 

• Update the nucleon (hadron) list and NN (hh) collision time list after current NN (hh) collision 
executed. 

• A new NN (hh) collision pair with least collision time is selected from the NN (hh) collision 
time list. 

• Repeat above six items until the NN (hh) collision time list is empty (hadronic freeze-out). 

3. A hadronic final state is eventually obtained for an A+B collision after the execution of parini_20b.f. 

4. usu.dat is an example for input file. The left panel of Fig. [1] shows the structure of PACIAE 2.0b 
program for the dynamical simulation of a nucleus-nucleus collision. 



C. PACIAE 2.0c structure 



PACIAE 2.0c program consists of paciae_20c.f, parini_20c.f, parcas_20c.f, sfm_20c.f, colaes_20c.f, had- 
cas_20c.f44, and p20c.f: 

1. paciae_20c.f plays double roles of the user program (an example) and the administrating an A+B 
collision event generation. 

2. parini_20c.f performs the parton initiation for an A+B collision. It is indeed a nucleon cascade 
process: 

• Initiation in the position and momentum spaces for an A+B collision is the same as that in 
PACIAE 2.0b. 

• The initial NN collision time list is also constructed as the same as that in PACIAE 2.0b. 

• A NN collision pair with least collision time is selected from the NN collision time list. 
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FIG. 4: (Color online) Charged particle transverse momentum distribution in 0-5% most central Pb+Pb collisions 
at ,/i^ r iv=2.76 TeV measured by ALICE and compared with the PACIAE 2.0b results. 



• This NN collision, if energy is large enough, executes by PYTHIA with the string fragmentation 
switched-off and the diquarks broken-up. The generated partons are filled in the parton list. 
If the energy is not large enough the current NN collision is dealt with usual two-body elastic 
collision [19| . 

• Update the nucleon list and NN collision time list. 

• A new NN collision pair with least collision time is selected from the NN collision time list. 

• Repeat above three items until the NN collision time list is empty (hadronic freeze-out). Then 
one obtains a partonic initial state (parton list) for an A+B collision. 

3. parcas_20c.f performs the parton rescattering for an A+B collision. 

4. sfm_20c.f executes the hadronization by the LUND string fragmentation regime (by calling p20c.f) 
for an A+B collision if switch parameter adjl(12)=0. 

5. coales_20c.f performs the hadronization for an A+B collision by the Monte Carlo coalescence model 
if switch parameter adjl(12)=l. 

6. hadcas_20c.f executes the hadron rescattering for an A+B collision. 

7. p20c.f is different from PYTHIA 6.4 in the addition of the mechanism for reduction of the strange 
quark suppression and of the requirement for charge conservation in the fragmentation of a string. 

The right panel of Fig. [T] shows the structure of PACIAE 2.0c program for the dynamical simulation 
of a nucleus-nucleus collision. Comparing left panel with right panel in Fig. [TJ one knows that in the 
PACIAE 2.0b model the parton initiation, parton rescattering, hadronization, and hadron rescattering 
arc performed for each hh collision pair independently. This is similar to the Monte Carlo Glauber model 
[3~0 ] despite that it is only a hadronic cascade model. PACIAE 2.0b characterizes the correlations among 
above processes for each hh collision pair, but there are no correlations among different hh collision pairs. 
This also means that in PACIAE 2.0b the nucleus-nucleus collision is dealt as a superposition of the 
nucleon-nucleon collisions in a sense. Hence the PACIAE 2.0b model presents the locality. Oppositely, 
the PACIAE 2.0c model characterizes the correlations among the different hh collision pairs in each 
process of the parton initiation, parton rescattering, hadronization, and hadron rescattering, but there 
are no correlations among above processes for each hh collision pair. PACIAE 2.0b and 2.0c are similar 
in the physical contents but are different in the topological structure. 

In order to decrease the line number in common block 'pyjets', a special common block 'sgam' is 
introduced. The photon is removed from 'pyjets' to 'sgam' after generation and is given specific flavor 
code as follows: 

• 22: hardonic decay photon, 

• 44: prompt direct photon, 

• 55: photon from parton-parton rescattering 

q + g -> q + i, Q + q-> 9 + 7, (82) 
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FIG. 5: (Color online) (a) charged particle pseudorapidity distribution in Au+Au collisions at - v /sjviv=200 GeV 
calculated by the PACIAE 2.0b (circles), PACIAE 2.0c (triangles), and PACIAE 2.0c_c (stars), (b) transverse 
momentum distribution. 



• 66: hardonic direct photon 

7T + 7T— >p + 7, 7T + /J ^ 7T + 7. (83) 

In default versions of PACIAE 2.0, the inelastic processes, qq -+ gg (denoted as process 6 in the 
program) and gg — > qq (process 7), are switched-off. If one wishes to include the inelastic process one 
has to comment out the statement of 'fsqq=0' and/or 'fsgg_l=0' in parcas_20a.f (or parcas_20b.f, or 
parcas_20c.f). 



IV. EXAMPLES CALCULATED 



The model parameters in PYTHIA 6.4 and PACIAE 2.0 were given default values according to the 
experimental measurement and/or the physical argument. However, in the specific calculations, a few 
sensitive parameters, as least as possible, should be tuned to a datum of the global measurable, such as 
the charged multiplicity or the charged particle rapidity density at mid-rapidity. The fitted parameters 
were then used to investigate the other physical observables. 

In the PACIAE 2.0a calculations the K factor was tuned to the charged particle pseudorapidity density 
at mid-pseudorapidity (|^| < 1) in the INEL (inelastic) pp collisions having at least one charged particle 
in the same region (INEL> O^i^i ) at y / s=0.9, 2.36, and 7 TeV [32[ . The results were shown in Tab. U 
together with the ALICE data [32J. This fitted K=1.9 (D=1.5, D means default value) was employed to 
calculate the charged particle pseudorapidity distribution and transverse momentum distribution (\r]\ < 
0.8) in the INEL pp collisions at -^=0. 9 TeV. These results were compared with the ALICE data |33l.l34| 
in Fig [5] One sees in this figure that the ALICE data are reproduced reasonably. 

We have used PACIAE 2.0b to calculate the 0-6% most central Au+Au collisions at ^sn jv=200 GeV 
and 0-5% most central Pb+Pb collisions at ^/sjvAf =2.76 TeV. In these calculations the K factor, the 
parameter f3 in LUND string fragmentation function, and the time accuracy At (the least time interval of 
two distinguishably consecutive collisions in the parton initiation stage) were tuned to the PHOBOS data 
of charged multiplicity Q and the ALICE data of charged particle pseudorapidity density at mid-rapidity 
(\rj\ < 0.5) [35|, respectively. The results were given in Tab. [II] The fitted parameters of K=1.7 (D=1.5), 
£=1.5 (D=0.58), and Ai=0.0001 for Au+Au as well as K=1.7, (3=1.5, and At =0.00004 for Pb+Pb were 
then employed to calculate the charged particle pseudorapidity distribution and transverse momentum 
distribution in Au+Au collisions and the transverse momentum distribution in Pb+Pb collisions by 
PACIAE 2.0b, respectively. These results were compared with the PHOBOS [IHI and the ALICE [33] 
data in Fig. [3] and |H respectively. We see in these figures that the experimental data are well described. 
The large fluctuation of the theoretical data at the tail of transverse momentum distribution in Fig. @] is 
because the total events of 1000 calculated is not enough. 

Similar calculations were also performed by PAICIA 2.0c and PAICIA 2.0c_c (which is the same as 
PACIAE 2.0c but hadronized by the Mote Carlo coalescence instead of the string fragmentation ). In 
the PAICIA 2.0c calculations K=1.7, (3—0.58, and A<=0.0001 were assumed for Au+Au collisions and 
K=1.7, /3=0.58, and A<=0.00001 for Pb+Pb collisions. The parameters of K=1.7 and Ai=0.0001 for 
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FIG. 6: (Color online) (a) Energy dependence of the participant scaled charged particle pseudorapidity density 
dN c h/drj/(N P art/'2) at midrapidity in 0-6% most central Au+Au collisions at RHIC energies and 0-5% Pb+Pb 
collisions at LHC energy, (b) Centrality dependence, where N nor indicates dN c h/ dr)/ (N par t/2) at largest N par t- 
In this figure the curves are fitted results and the open symbols are results from PACIAE 2.0b calculations. 



Au+Au collisions as well as K=1.7 and At—0. 000055 for Pb+Pb collisions were assumed in PAICIA 
2.0c_c calculations. The results of charged particle multiplicities were compared with the results from 
PACIAE 2.0b calculations and the experimental data in Tab. |TTJ Fig. [5] gave the comparison of the 
charged particle pseudorapidity (panel (a)) and transverse momentum (panel (b)) distributions with the 
corresponding PACIAE 2.0b results above. We see in this figure that PAICIA 2.0c is also able to describe 
the experimental data. However, the discrepancy between PAICIA 2.0c and PAICIA 2.0c_c is visible. 

Recently, PHOBOS investigated the energy dependence of the participant scaled charged particle 
pseudorapidity density at mid-rapidity dN c f l /dr)/(N par t/2) in the relativistic pp and A+B collisions. They 
obtained the function of dN c h/ drj/ (N part /2) vs. ^/snn by fitting the experimental data spanned from 10 
GeV to 7 TeV for pp and from 2 GeV to 200 GeV for A+B collisions [lj (Before the submission of this 
paper we knew that PHOBOS just published their paper in Phys. Rev. C where they added a fitting of the 
experimental data from 19.6 GeV to 2.76 TeV with a power law in ^/snn for A+B collisions.). Assuming 
further that the dependence of dN ca /dri/(N par t/2) on energy and centrality (denoted by Glauber model 
Np ar t) can be factorized, then they obtained the function of dN c h/ drj/{N part /2) vs. ^smn and the 
function of dN ch /dr)/(N part /2) vs. N part . 

We used the geometric number of participant nucleons instead of the Glauber model N part to fit first 
the PHOBOS data of 0-6% most central Au+Au collisions at RHIC energies [2 and ALICE datum of 0- 
5% most central Pb+Pb collisions at v /i^T=2.76 TeV [35] (cf. Fig. [6] (a)). Then we fitted the PHOBOS 
centrality dependence data [l[ with normalization to the datum of highest centrality (i. e. N nor in the 
figure) and ALICE centrality dependence data [35[ with similar normalization (cf. Fig.HJb)). The fitted 
functions were obtained 

z = f{x)-^- v (84) 

9 '{Umax ) 

f(x) = 0.2 In(a;) + 0.001788(ln(x)) 2 + 0.01244(ln(a;)) 3 + 0.8916, (85) 
g(y) = -0.2217 + 0.3088y (1/3) - 0.01906j/ 2/3) , (86) 
where z,x, and y denoted dN c h/ 'dr]/(N par t/2), tJsnn, and Npart, respectively. 

a 



We fitted the PHOBOS pp data jlj] and ALICE pp data [32j of dN cn /dr] by the function as the same 
as Eq. (JHU). The new four coefficients of 1.725, -0.1312, 0.007777, and -4.3842 were obtained (cf. Fig. [7J. 

With above fitted relations one may be able to first predict dN c h/dr]/(N pa rt/2) for un-measured pp 
and/or A+B collisions in the RHIC and LHC energy region. Then tuning the sensitive model parameters 
(such as K factor, j3, and/or At above) to this dN ca / dr\ / (N part / '2) one may also be able to predict other 
observables such as pseudorapidity distribution, transverse momentum distribution, etc.. 
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FIG. 7: (Color online) Energy dependence of the participant scaled charged particle pseudorapidity density 
dN c h/dri at midrapidity in pp collisions at RHIC and LHC energies. The curve is fitted results and the open 
symbols are results from PACIAE 2.0a calculations. 



APPENDIX: PACIAE 2.0 USERS' GUIDE 



1. Selection of the processes 

A switch variable (nchan) is introduced as follows: 



• nchan=0: 


Inelastic (INEL), 


• nchan=l: 


Non Single Diffractive (NSD), 


• nchan=2: 


qq — > 7*/Z° used to generate Drell-Yan process, 


• nchan=3: 


3 ftp production, 


• nchan=4: 


heavy-flavor production, 


• nchan=5: 


direct photon. 


• nchan=6: 


soft processes only, 


• nchan=7: 


default PYTHIA. 


2. Run PACIAE 



• Compile each of the FORTRAN program packets in 20a. tar. gz or 20b. tar. gz or 20c. tar. gz and link 
them together forming an executable file. 

• Set the parameters in input file of usux.dat for PACIAE 2.0a or usu.dat for PACIAE 2.0b and 
PACIAE 2.0c properly. 

• Run the executable file. 
3. Parameters 

The parameters introduced in PACIAE 2.0, except those in PYTHIA 6.4, can be found in the read 
statements in user program of paciae_20a.f (paciae_20b.f, paciae_20c.f). Meaning of parameter can be 
found in the program packets via searching its name. Those parameters not described in the program 
packets and the array adjl(40) are explained as follows: 
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• nout: internal output per 'nout' events. 

• psno=l: systematic sampling method for impact parameter in A+B collisions, 
psno=2: randomly sampling method for impact parameter. 

To run an event with a given impact parameter, set 'bmin'='bmax' in usu.dat. 

• adjl(i), i= 

1: K factor in parton cascade. 

2: a s , effective coupling constant in the parton rescattering (cascade). 

3: parameter 'tcut' (i cut ) introduced in the integration of parton-parton differential 

cross section in the subroutine 'fsig' in parcas_20a.f (parcas_20b.f, parcas_20c.f). 
4: parameter 'idw', the number of intervals in the numerical integration 

in parcas_20a.f (parcas_20b.f, parcas_20c.f). 
5: =1 with nuclear shadowing, 

=0 without nuclear shadowing. 
6: parameter a (parj(41) in PYTHIA 6.4) in the LUND string fragmentation function. 
7: parameter j3 (parj(42) in PYTHIA 6.4) in the LUND string fragmentation function. 
8: mstp(82) in PYTHIA 6.4 

=1, with hard interactions, 

=0: without hard interactions (simple two-string model, soft only). 
9: parp(81) in PYTHIA 6.4 (default=1.4 GeV/c), effective minimum transverse momentum for 

multiple interactions when mstp(82)=l. 
10: K factor (parp(31) in PYTHIA 6.4). 

11: time accuracy used in the hadron cascade (hadcas_20a.f, hadcas_20b.f, and hadcas_20c.f). 
12: model for hadronization 

=0 string fragmentation, 

=1 Monte Carlo coalescence model. 
13: dimension of meson table considered. 
14: dimension of baryon table considered. 
15: string tension. 

16: number of loops in the deexcitation of energetic quark in the Monte Carlo coalescence model. 
17: the threshold energy in the deexcitation of energetic quark in the Monte Carlo coalescence model. 
18: =0 without Pauli blocking in the parton cascade, 

=1 with Pauli blocking. 
19: time accuracy used in the parton cascade (parcas_20a.f, parcas _20b.f, and parcas_20c.f). 
20: =0 using LO pQCD parton-parton cross section in the parton rescattering (cascade), 

=1 using keeping only leading divergent terms in the LO pQCD parton-parton cross section, 

=2 the same as but flat scattering angle distribution is assumed, 

=3 the same as 1 but flat scattering angle distribution is assumed. 
21: =0 without phase space constraint in the Monte Carlo coalescence model, 

=1 with phase space constraint. 
22: critical value of the product of radii in the position and momentum phase spaces (4 is assumed). 
23: =0 LUND string fragmentation function is used in the subroutine for deexcitation of the 

energetic quark in the Monte Carlo coalescence model, 

=1 Field- Feynman parametrization function is used. 
24: the virtuality cut ('tlO') in the time-like radiation in the parton rescattering (cascade). 
25: Aqcd m the parton cascade. 
26: number of random number thrown away. 
27: largest momentum allowed for particle. 

28: concerned to the largest position allowed for particle, see program for the detail. 

29: width of two dimension Gaussian distribution sampling p x and p y of the produced quark 

pair in the deexcitation of energetic quark in the Monte Carlo coalescence model. 
30: maximum -p\ in above two dimension Gaussian distribution. 
31: parj(l) in PYTHIA 6.4. 
32: parj(2) in PYTHIA 6.4. 
33: parj(3) in PYTHIA 6.4. 
34: parj(21) in PYTHIA 6.4. 

35: mstp(91) in PYTHIA 6.4, parton transverse momentum (k±) distribution inside hadron; 
=1, Gaussian; 
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=2, exponential 

36: =0 without phenomenological parton energy loss in the parton rescattering (cascade), 

=1 with phenomenological parton energy loss. 
37: the coefficient in phenomenological parton energy loss. 
38: pt cut in phenomenological parton energy loss. 
39: width of Gaussian k± distribution in hadron if mstp(91)=l, 

width of exponential kj_ distribution in hadron if mstp(91)=2. 
40: decide where the run stopped: 

=1, after parton initiation; 

=2, after parton rescattering; 

=4, after hadron rescattering. 

4. Output files 

• rmsO.out is an output for the input parameters. 

• rms.out is the main output hie. 

• main. out is the PYTHIA standard output. 

• oscar.out is the OSCAR standard output if nosc=2 or 3. 

• Encourage user to write his/her own output hie and user program. 

5. Postscript 

• The PACIAE 2.0 program is free for the person who is interested. However, it is protected by GPL, 
no merchandise use please. 

• Welcome reporting any bugs to yanyl@ciae.ac.cn, zhoudm@phy.ccnu.edu.cn , and/or 
sabh@ciae.ac.cn . 
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